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Abstract 

We calculate LSP-nucleus elastic scattering cross sections using some representative in¬ 
put in the restricted SUSY parameter space. The coherent matrix elements are computed 
throughout the periodic table while the spin matrix elements for the proposed target 

which has a rather simple nuclear structure. The results are compared to those given from 
other cold dark matter detection targets. 
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1 Introduction 


There are many arguments supporting the fact that, the cold dark matter of the universe, i.e. 
its component which is composed of particles which were non-relativistic at the time of structure 
formation, is at least 60%. There are two interesting cold dark matter candidates: i) Massive 
Compact Halo Objects (MACHO’s) and ii) Weak Interacting Massive Particles (WIMP’s). In the 
present work we discuss a special WIMP candidate connected with the supersymmetry, i.e. the 
lightest supersymmetric particle (LSP). 

We examine the possibility to directly detect the LSP via the recoiling of a nucleus (A,Z) in the 
elastic scattering process: y + (A, Z) ^ x + Z)* (y denotes the LSP). In this investigation, 

we proceed with the following steps: 

1) We write down the effective Lagrangian at the elementary particle (quark) level obtained in 
the framework of supersymmetry as described in refs. P, ^ 

2) We go from the quark to the nucleon level using an appropriate quark model for the nucleon. 
Special attention in this step is paid to the scalar couplings, which dominate the coherent part of 
the cross section and the isoscalar axial currents which depend on the assumed quark model. 

3) We compute the relevant nuclear matrix elements using as reliable as possible many boby 
nuclear wave function hopping that, by putting as accurate nuclear physics input as possible, we 
will be able to constrain the SUSY parameters as much as possible. 

4) We calculate the modulation of the cross sections due to the earth’s revolution around the sun 
and the motion of the sun around the center of the galactic disc by a folding procedure. 

There are many popular targets for LSP detection as 
etc. Among them Pb has been recently proposed as an important detector, since its spin 
matrix element, especially the isoscalar one, does not exhibit large quenching as that of the light 
and up to now much studied and nuclei. 

Our purpose is to calculate LSP-nucleus scattering cross section using some representative 
input in the restricted SUSY parameter space, to compute the coherent LSP-nucleus scattering 
cross sections throughout the periodic table and study the spin matrix elements of since this 

target, in addition to its experimental qualifications, has the advantage of a rather simple nuclear 
structure. We compare our results to those given from other proposed cold dark matter detection 
targets |0] and present results calculated by using new imput SUSY parameters ^ obtained in a 
phenomenologically allowed parameter space. 


2 Momentum transfer dependence of the total cross sec¬ 
tions 

The total cross section of the LSP-nucleus elastic scattering can be written as 0] 
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where mi, rUp is the LSP, proton masses, r] = mi/nipA and cxo — 0.77x 10“^®cm^. The parameters 
/s) fvy /a) with r = 0,1 an isospin index, describe the scalar, vector and axial vector couplings 
and are determined by a specihc SUSY mechanism. H The momentum transfer enters via u as 
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(b is the harmonic oscillator parameter). In eq. (|l|), the terms in the square brackets describe the 
coherent cross section and involve the integrals Ip{u) {p = 0,1) while the other terms describe the 
spin dependence of the cross section by means of the integrals Ippiu) (p, p' = 0,1, isospin indices) 
which are normalized so as /p(0) = 1 and /pp'(O) = 1. For their dehnitions see ref. ^ where they 
are calculated by using realistic nuclear form factors. 

In ref. [|] we studied the variation of Ip as functions of the LSP mass. We found that u can 
be quite big for large mass of the LSP and heavy nuclei even though the energy transfer is small 
(< lOOiPeU). The total cross section can in such instances be reduced by a factor of about hve. 

The spin matrix element of heavy nuclei like has not been evaluated, since one expects the 
relative importance of the spin versus the coherent mode to be more pronounced on light nuclei. 
However, the spin matrix element in the light isotopes is quenched, while that of does 

not show large quenching. For this feature, we recently proposed nucleus as an important 

candidate in the LSP detection. ^ Furthermore, this nucleus has some additional advantages as: 
i) it is believed to have simple structure (one 2pl/2 neutron hole outside the doubly magic nucleus 
208 p^) g^j^d ii) it has low angular momentum and therefore only two multipoles A = 0 and A = 2 
with a J-rank of k = 1 can contribute even at large momentum transfers. 

The dependence of the total cross section on the momentum transfer for has been studied 
in ref. Q where the variation of the integrals Ipp/{u), which describe the spin part of ‘^^'^Pb, has 
also been investigated. 


3 Folding of the cross section with a Maxwellian distribu¬ 
tion 


Due to the revolution of the earth around the sun, the LSP-nucleus cross sections appear modihed. 
This effect is studied by folding the total cross section eq. (|^) with an appropriate distribution. If 
we assume a Maxwell type distribution, which is consistent with the velocity distribution of LSP 
into the galactic halo, the counting rate for a target with mass m takes the form [Q 
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with p(0) = 0.3GeU/cm^, the LSP density in our vicinity, and < S > given by 
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Table 1: Comparison of the static spin matrix elements for three typical nuclei. 


Component 

207 pu 

''^Gegi2+ 

DZi/2+ 

“B?(0) 

0.231 

1.005 

0.204 

Bi(0)r!o(0) 

-0.266 

-1.078 

-0.202 

B§(0) 

0.305 

1.157 

0.201 


< S > = 


The parameters Jo, Jp, Jpa describe the scalar, vector and spin part of the velocity averaged 
counting rate, respectively. They are functions oi (3 q = vq/ c, 5 = 2vi/vq and mq, where uq is the 
velocity of the sun around the galaxy, vi the velocity of the Earth around the sun and uq is given 
by an expression like that of eq. (|^) with f3 —*• Pq. Since 5 ~ 0.27 << 1, we can expand the 
J-integrals in powers of 5 and retain terms up to linear in 5. Thus, for each mechanism (vector, 
scalar, spin) we obtain two integrals associated with I = 0 and / = 1. The most important 
integrals are studied in ref. By exploiting the above expansion of K-integrals, the counting rate 
can be written in the form 
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where a is the phase of the earth’s orbital motion and {dN/dt)^ is the rate obtained from the 
I = 0 multipole, h represent the amplitude of the oscillation, i.e. the ratio of the component of 
the multipole / = 1 to that of / = 0. 


4 Results and discussion 


The three basic ingredients of our calculation were the input SUSY parameters, a quark model 
for the nucleon and the structure of the nuclei involved. The input SUSY parameters used for 
the results presented in Tables 1 and 2 have been calculated in a phenomenologically allowed 
parameter space (cases #1, #2, ^3) as explained in ref. [Q 

For the coherent part (scalar and vector) we used realistic nuclear form factors and studied 
three nuclei, representaves of the light, medium and heavy nuclear isotopes {Ca, Ge and Pb). In 
Table 1 we show the results obtained for three different quark models denoted by A (only quarks 
u and d) and B, C (heavy quarks in the nucleon). We see that the results vary substantially and 
are sensitive to the presence of quarks other than u and d into the nucleon. 

The spin contribution, arising from the axial current, was computed in the case of system. 
For the isovector axial coupling the transition from the quark to the nucleon level is trivial (a factor 











Table 2: The quantity < dN/dt >o in y ^Kg ^ and the oscillator parameter h for the coherent 
vector and scalar contributions. 



Vector Contribution 

Scalar Contribution 


< dN/dt >0 

h 


< dN/dt >0 


h 

Case 

(xlO-Q 

Model A 

Model B 

Model C 



#1 

0.264 

0.029 

0.151 X 10-^ 

0.220 

0.450 

-0.002 

Pb 

#2 

0.162 

0.039 

0.410 X 10-1 

142.860 

128.660 

0.026 


#3 

0.895 

0.038 

0.200 X 10-3 

0.377 

0.602 

-0.001 


#1 

0.151 

0.043 

0.779 X 10-^ 

0.120 

0.245 

0.017 

Ge 

#2 

0.053 

0.057 

0.146 X 10-1 

51.724 

46.580 

0.041 


#3 

0.481 

0.045 

0.101 X 10-3 

0.198 

0.316 

0.020 


#1 

0.079 

0.053 

0.340 X 10-1 

0.055 

0.114 

0.037 

Ca 

#2 

0.264 

0.060 

0.612 X 10-2 

22.271 

20.056 

0.048 


#3 

0.241 

0.053 

0.435 X 10-1 

0.090 

0.144 

0.038 


Table 3: The spin contribution in the LSP Pb scattering for two cases: EMC data and NQM 
Model. The LSP mass is mi = 126, 27, 102 GeV for #1, #2, #3 respectively. 



EMC DATA 


NQM 

MODEL 

Solution 

< dN/dt >0 (y-iATg-i) 

h 

< dN/dt >0 {y '^Kg i) h 

#1 

0.285 X 10-2 

0.014 

0.137 X 10-2 

0.015 

#2 

0.041 

0.046 

0.384 X 10-2 

0.056 

#3 

0.012 

0.016 

0.764 X 10-2 

0.017 


of qa = 1.25). For the isoscalar axial current we considered two possibilities depending on the 
portion of the nucleon spin which is attributed to the quarks, indicated by EMC and NQM. [Q 
The ground state wave function of obtained by diagonalizing the nuclear Hamiltonian 

in a 2h-lp space which is standard for this doubly magic nucleus. The momentum dependence of 
the matrix elements was taken into account and all relevant multipoles were retained (here only 
A = 0 and A = 2). 

In Table 3, we compare the spin matrix elements at g = 0 for the most popular targets 
considered for LSP detection and We see that, even though the spin matrix 

elements 12^ are even a factor of three smaller than those for obtained in ref. [0 (see Table 
1), their contribution to the total cross section is almost the same (see Table 4) for LSP masses 
around 100 GeH. Our hnal results for the quark models (A, B, C, NQM, EMC) are presented in 
Tables 3 and 4 for SUSY models #l-#3 P] and Tables 4 and 5 for SUSY models ||| 


5 Conclusions 


In the present study we found that for heavy LSP and heavy nuclei the results are sensitive to 
the momentum transfer as well as to the LSP mass and other SUSY parameters. From the Tables 
3 and 4 we see that, the results are also sensitive to the quark structure of the nucleon in the 
following sense: i) The coherent scalar (associated with Higgs exchange) for model A (u and d 
quarks only) is comparable to the vector coherent contribution. Both are at present undetectable. 
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Table 4: Ratio of spin contribution {^^'^Pb/'^^Ge) at the relevant momentum transfer with the 
kinematical factor 1/(1 + rff, rj = mi/Arrip. 


Solution 

#1 

#2 

#3 

#4 

#5 

#6 

#7 

#8 

#9 

mi (GeV) 

126 

27 

102 

80 

124 

58 

34 

35 

50 

NQM 

0.834 

0.335 

0.589 

0.394 

0.537 

0.365 

0.346 

0.337 

0.417 

EMC 

0.645 

0.345 

0.602 

0.499 

0.602 

0.263 

0.341 

0.383 

0.479 


Table 5: The same parameters as in Tables 1 and 2. The LSP mass is mi = 80, 124, 58, 34, 35, 
50 GeV for the cases ^(^4 — ^9 respectively. Cases 7(^8, #9 are no-scale models. The values of 
< dN/dt >0 for Model A and the Vector part must be multiplied by xl0“^. 




Scalar 

Part 


Vector 

Part 

Spin Part 



\ dt /() 


h 

\ d.t / C\ 

h 


\ dt /n 


h 


A 

B 

C 



EMC 


NQM 



#4 

0.03 

22.9 

8.5 

0.003 

0.04 

0.054 

0.8010- 


0.1610' 


0.015 

#5 

0.46 

1.8 

1.4 

-0.003 

0.03 

0.053 

0.3710- 

-3 

0.9110' 

-3 

0.014 

#6 

0.16 

5.7 

4.8 

0.007 

0.11 

0.057 

0.4410- 

-3 

0.1110' 

-2 

0.033 

#7 

4.30 

110.0 

135.0 

0.020 

0.94 

0.065 

0.67 


0.87 


0.055 

#8 

2.90 

73.1 

79.8 

0.020 

0.40 

0.065 

0.22 


0.35 


0.055 

#9 

2.90 

1.6 

1.7 

0.009 

0.95 

0.059 

0.29 


0.37 


0.035 


ii) For models B and C (heavy quarks in the nucleon) the coherent scalar contribution is dominant. 
Detectable rates < dN/dt >o> 100 y~^Kg~^ are possible in a number of models with light LSP. 

The spin contribution is sensitive to the nuclear structure. It is undetectable if the LSP is 
primarily a gaugino. The folding of the event rate with the velocity distribution provides the 
modulation effect h. In all cases it is small, less than ±5%. 
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